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Abstract 

We present a systematic method for reducing an arbitrary one— loop TV— point massless 
Feynman integral with generic 4— dimensional momenta to a set comprised of eight fundamen- 
tal scalar integrals: six box integrals in D — 6, a triangle integral in D = 4, and a general 
two— point integral in D space time dimensions. All the divergences present in the original inte- 
gral are contained in the general two— point integral and associated coefficients. The problem of 
vanishing of the kinematic determinants has been solved in an elegant and transparent manner. 
Being derived with no restrictions regarding the external momenta, the method is completely 
general and applicable for arbitrary kinematics. In particular, it applies to the integrals in 
which the set of external momenta contains subsets comprised of two or more collinear mo- 
menta, which are unavoidable when calculating one— loop contributions to the hard— scattering 
amplitude for exclusive hadronic processes at large momentum transfer in PQCD. The iterative 
structure makes it easy to implement the formalism in an algebraic computer program. 
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1 Introduction 



Scattering processes have played a crucial role in establishing the fundamental interactions of nature. 
They represent the most important source of information on short-distance physics. With increasing 
energy, multiparticle events are becoming more and more dominant. Thus, in testing various aspects 
of QCD, the high-energy scattering processes, both exclusive and inclusive, in which the total number 
of particles (partons) in the initial and final states is N > 5, have recently become increasingly 
important. 

Owing to the well-known fact that the LO predictions in perturbative QCD (PQCD) do not have 
much predictive power, the inclusion of higher-order corrections is essential for many reasons. In 
general, higher-order corrections have a stabilizing effect, reducing the dependence of the LO predic- 
tions on the renormalization and factorization scales and the renormalization scheme. Therefore, to 
achieve a complete confrontation between theoretical predictions and experimental data, it is very 
important to know the size of radiative corrections to the LO predictions. 

Obtaining radiative corrections requires the evaluation of one-loop integrals arising from the 
Feynman diagram approach. With the increasing complexity of the process under consideration, 
the calculation of radiative corrections becomes more and more tedious. Therefore, it is extremely 
useful to have an algorithmic procedure for these calculations, which is computerizable and leads to 
results which can be easily and safely evaluated numerically 

The case of Feynman integrals with massless internal lines is of special interest, because one 
often deals with either really massless particles (gluons) or particles whose masses can be neglected in 
high— energy processes (quarks). Owing to the fact that these integrals contain IR divergences (both 
soft and collinear), they need to be evaluated in an arbitrary number of space-time dimensions. As it 
is well known, in calculating Feynman diagrams mainly three difficulties arise: tensor decomposition 
of integrals, reduction of scalar integrals to several basic scalar integrals and the evaluation of a set 
of basic scalar integrals. 

Considerable progress has recently been made in developing efficient approaches for calculating 
one— loop Feynman integrals with a large number (N > 5) of external lines [Tl 01151111 151 Ifil 171151 1§1 1 lOj. 
Various approaches have been proposed for reducing the dimensionally regulated (N > 5)— point 
tensor integrals to a linear combination of iV— and lower— point scalar integrals multiplied by tensor 
structures made from the metric tensor g^ v and external momenta P3 Gl El El di| ■ It has also 
been shown that the general (N > 5)— point scalar one— loop integral can recursively be represented 
as a linear combination of (N — 1)— point integrals provided the external momenta are kept in 
four dimensions El El El El El • Consequently, all scalar integrals occurring in the computation 
of an arbitrary one— loop (N > 5)— point integral can be reduced to a sum over a set of basic 
scalar box (N = 4) integrals with rational coefficients depending on the external momenta and 
the dimensionality of space— time. Despite, the considerable progress, the developed methods still 
cannot be applied to all cases of practical interest. The problem is related to vanishing of various 
relevant kinematic determinants. 

As far as the calculation of one-loop (N > 5)-point massless integrals is concerned, the most 
complete and systematic method is presented in |7J. It does not, however, apply to all cases of 
practical interest. Namely, being obtained for the non-exceptional external momenta it cannot be, 
for example, applied to the integrals in which the set of external momenta contains subsets comprised 
of two or three collinear on-shell momenta. The integrals of this type arise when performing the 
leading-twist NLO analysis of hadronic exclusive processes at large— momentum transfer in PQCD. 

With no restrictions regarding the external kinematics, in this paper we formulate an efficient, 
systematic and completely general method for reducing an arbitrary one— loop N— point massless 
integral to a set of basic integrals. Although the method is presented for massless case, the gener- 
alization on massive case is straightforward. The main difference between the massive and massless 
cases manifests itself in the basic set of integrals, which in former case is far more complex. Among 
the one— loop Faynman integrals there exist both massive and massless integrals for which the ex- 
isting reduction methods break down. The massless integrals belonging to this category are of more 
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practical interest at the moment, so in this paper we concentrate on massless case. 

The paper is organized as follows. Section 2 is devoted to introducing notation and to some 
preliminary considerations. In Sec. 3, for the sake of completenes, we briefly review a tensor 
decomposition method for N— point tensor integrals which was originally obtained in Ref. [2]. In 
Sec. 4 we present a procedure for reducing one— loop N— point massless scalar integrals with generic 
4— dimensional external momenta to a fundamental set of integrals. Since the method is closely 
related to the one given in E], similarities and differences between the two are pointed out. 
Being derived with no restrictions to the external momenta, the method is completely general and 
applicable for arbitrary kinematics. Section 5 contains considerations regarding the fundamental set 
of integrals which is comprised of eight integrals. Section 6 is devoted to some concluding remarks. 
In the Appendix A we give explicit expressions for the relevant basic massless box integrals in D = 6 
space— time dimensions. These integrals constitute a subset of the fundamental set of scalar integrals. 
As an illustration of the tensor decomposition and scalar reduction methods, in the Appendix B we 
evaluate an one— loop 6— point Feynman diagram with massless internal lines, contributing to the 
NLO hard— scattering amplitude for 77-^ tt + tt~ exclusive reaction at large momentum transfer in 
PQCD. 



2 Definitions and general properties 

In order to obtain one— loop radiative coorections to physical processes in massless gauge theory, 
the integrals of the following type are required: 

i N CD- -To -V) = (u 2 ) 2 - D/2 ( d ° l lt±1 ' " ■ m 

V^^'iWi-^i J (2tt) d A 1 A 2 ---A n W 

This is a rank P tensor one— loop iV— point Feynman integral with massless internal lines in D— di- 
mensional space— time, where pi, (i — 1, 2, N) are the external momenta, I is the loop momentum, 
and /i is the usual dimensional regularization scale. 

The Feynman diagram with N external lines, which corresponds to the above integral, is shown 
in Fig. ^ For the momentum assignments as shown, i.e. with all external momenta taken to be 
incoming, the massless propagators have the form 

Ai = (l + ri ) 2 +ie i = 1, - ■ ■ , N , (2) 

where the momenta r, are given by r% = pi + r^_i for i from 1 to N, and ro = rjy. The quantity 
ie(e > 0) represents an infinitesimal imaginary part, it ensures causality and after the integration 
determines the correct sign of the imaginary part of the logarithms and dilogarithms. It is customary 
to choose the loop momentum in a such a way that one of the momenta vanishes. However, for 
general considerations, such a choice is not convenient, since by doing so, one loses the useful 
symmetry of the integral with respect to the indices 1, • • • , N. 
The corresponding scalar integral is 

iSiD-M^^f-^j^^^-. (3) 

If P + D — 2N > 0, the integral JIJ is UV divergent. In addition to UV divergence, the integral can 
contain IR divergence. There are two types of IR divergence: collinear and soft. A Feynman diagram 
with massless particles contains a soft singularity if it contains an internal gluon line attached to two 
external quark lines which are on mass— shell. On the other hand, a diagram contains a collinear 
singularity if it contains an internal gluon line attached to an external quark line which is on mass 
—shell. Therefore, a diagram containing a soft singularity at the same time contains two collinear 
singularities, i.e. soft and collinear singularities overlap. 
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Figure 1: One-loop N-point diagram 



When evaluating Feynman diagrams, one ought to regularize all divergences. Making use of the 
dimensional regularization method, one can simultaneously regularize UV and IR divergences, which 
makes the dimensional regularization method optimal for the case of massless field theories. 

The tensor integral Q is, as it is seen, invariant under the permutations of the propagators Ai, 
and is symmetric with respect to the Lorentz indices Lorentz covariance allows the decomposition 
of the tensor integral £Q) in the form of a linear decomposition consisting of the momenta pi and the 
metric tensor g^ v . 



3 Decomposition of tensor integrals 

Various approaches have been proposed for reducing the dimensionally regulated TV— point tensor 
integrals to a linear combination of TV— and lower— point scalar integrals multiplied by tensor struc- 
tures made from the metric tensor g^ v and external momenta. In this section we briefly review the 
derivation of the tensor reduction formula originally obtained in Ref. [2] • 

For the purpose of the following discussion, let us consider the tensor integral 

I N (D-iv-\) = (u 2 ) 2 - D/2 [ hulll l m a") 

Vi-m-pv^' — > J (2tt) d A^Az 2 ■ ■ ■ A v ^ ' 

and the corresponding scalar integral 

gmwi-o^/g, ^ 1 ,,,,., • w 

The above integrals represent generalizations of the integrals and ©, in that they contain 
arbitrary powers V{ £ N of the propagators in the integrand, where {vi} is the shorthand notation 
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for (yi, • • • , vn)- Also, for notational simplicity, the external momenta are omitted from the argument 
of the integral. 

The Feynman parameter representation of the tensor integral I^—ap (D> {^i})> given in fllf. which 
is valid for arbitrary values of N, P, and Vi(> 0), for the values of D for which the remaining 
integral is finite, and the T— function does not diverge is given by 



/l., P (AW) = 



(4tt) 2 

r(Ei^-J>/2-fc) 
2 fc in 4 r(^)] 



(4V) 



2\2-£>/2 



E {[5fh? 1 ---[r J vr} Ml 



2fc + Ej.;=.P 



(-1) 



fc+D/2-Eii/i 



(6) 



where • • • [?'jvF jv j>i---/i P represents a symmetric (with respect to //i • • • /ip) combination of 

tensors, each term of which is composed of k metric tensors and ji external momenta r^. Thus, for 
example, 

As for the integral representation of the corresponding scalar integral the result is of the form 



IS 1 (DM) 



(4tt) ; 



- {Ann 



2\2-D/2 



r(EL»i-D/2) 



N 



i<3 



(7) 



Now, on the basis of J7J, (JSJ can be written in the form 

<.., p (aw)= E WW 1 



[r N ] jN } 



(4V) 



2\P-k 



{-2f 



' N 

n 



i>0 



J^P + 2(P-fc);{^+ Jl }). 



(8) 



This is the desired decomposition, of the dimensionally regulated N— point rank P tensor integral. It 
is originally obtained in Ref. [2J. Based on JSJ), any dimensionally regulated N— point tensor integral 
can be expressed as a linear combination of N— point scalar integrals multiplied by tensor structures 
made from the metric tensor and external momenta. Therefore, with the decomposition JSJ, the 
problem of calculating the tensor integrals has been reduced to the calculation of the general scalar 
integrals. 

It should be pointed out that among the tensor reduction methods presented in the literature 
one can find methods, e.g. |7j which for N > 5 completely avoid the terms proportional to the 
metric tensor g^" . Compared with the method expressed by iJSJ, these reduction procedures lead 
to decomposition containing a smaller number of terms. The methods of this type are based on 
the assumption that for N > 5, one can find four linearly independent 4- vectors forming a basis 
of the 4-dimensional Minkowski space, in terms of which the metric tensor can then be expressed. 
This assumption is usually not realized when analyzing the exclusive processes at large- momentum- 
transfer (hard scattering picture) in PQCD. Thus, for example, in order to obtain the next-to-leading 
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order corrections to the hard-scattering amplitude for the proton-Compton scattering, one has to 
evaluate one-loop N — 8 diagrams. The set of external momenta contains two subsets comprised of 
three collinear momenta (representing the proton) . The kinematics of the process is thus limited to 
the 3-dimensional subspace. If this is the case, the best way of doing the tensor decomposition is 
the one based on formula (JSJ , regardless of the fact that for large N the number of terms obtained 
can be very large. 

As is well known, the direct evaluation of the general scalar integral (JHJ) (i.e. 10) represents a non- 
trivial problem. However, with the help of the recursion relations, the problem can be significantly 
simplified in the sense that the calculation of the original scalar integral can be reduced to the 
calculation of a certain number of simpler fundamental (basic) intgerals. 

4 Recursion relations for scalar integrals 

Recursion relations for scalar integrals have been known for some time 21 EH U\- However, 
as it turns out, the existing set of relations that can be found in the literature is not sufficient to 
perform the reduction procedure completely, i.e. for all one-loop integrals appearing in practice. 
The problem is related to vanishing of various kinematic determinants, it is manifest for the cases 
corresponding to N > 6, and it is especially acute when evaluating one-loop Feynman integrals 
appearing in the NLO analysis of large momentum transfer exclusive processes in PQCD. As is well 
known, these processes are generally described in terms of Feynman diagrams containing a large 
number of external massless lines. Thus, for example, for nucleon Compton scattering is N = 8. 
A large number of external lines implies a large number of diagrams to be considered, as well as 
a very large number of terms generated when performing the tensor decomposition using (JHJ. In 
view of the above, to treat the Feynman integrals (diagrams) with a large number of external lines 
the use of computers is unavoidable. This requires that the scalar reduction procedure be generally 
applicable. It is therefore absolutely clear that any ambiguity or uncertainty present in the scalar 
recursion relations constitutes a serious problem. The method presented below makes it possible 
to perform the reduction completely regardless of the kinematics of the process considered and the 
complexity of the structure of the contributing diagrams. 

For the reason of completeness and clearness of presentation and with the aim of comparison with 
the already existing results, we now briefly present a few main steps of the derivation of recursion 
relations. It should be pointed out that the derivation essentially represents a variation of the 
derivation originally given in yy. 

Recursion relations for scalar integrals are obtained with the help of the integration-by-parts 
method [SJEUE]- Owing to translational invariance, the dimensionally regulated integrals satisfy 
the following identity: 



where Zi (i = • • • N) are arbitary constants, while Ai are the propagators given by (0). The identity 
@ is a variation of the identity used in [S], where it was assumed that — 0. Performing the 
differentiation, expressing scalar products in the numerator in terms of propagators A i: choosing 
zq = Ylz=i z ii (which we assume in the following) and taking into account the scalar integral JSJ, 
the identity @ leads to the relation 




(9) 



N / N 





j=l \i=l 



N N 



]T ZiVjlq (D; {v k + S kj - 8 kl }) - (D - VjW${D; {v k }), 



(10) 




5 



where Sij is the Kronecker delta symbol. In arriving at i|l(Jf) . it has been understood that 

Iq{D; vi, ■ ■ - .fi-ijO, vt + i, ■■■,v N ) = I^ 1 (D; v u ■ ■ ■ , vi_ x , n + x, ■ ■ ■ ,v N ). (11) 

The relation (|1U|) represents the starting point for the derivation of the recursion relations for scalar 
integrals. 

We have obtained the fundamental set of recursion relations by choosing the arbitrary constants 
Zi so as to satisfy the following system of linear equations: 

E- M-r j fzi = C, j = l,...,N, (12) 

* '2 — 1 

where C is an arbitrary constant. Introducing the notation r%j = (r 4 — rj) 2 , the system (11211 may be 
written in matrix notation as 



t 


ri2 ■ 


■ r 1N 


\ 


/ 


Z\ 


\ 


( C \ 




• 


■ r 2 N 




Z2 


c 


\ r 1N 


T2N ■ 


■ 


/ 


V 


ZN 


J 


V c ) 



It should be pointed out that the expression of the type ifTHjl and the system of the type l(T3|l , for 
the case of massive propagators (Aj = (/ + r 4 ) 2 — m 2 + ie), see Ref. [S], can simply be obtained from 
the relation given above by making a change rij — » rij — m 2 — m 2 . Consequently, considerations 
performed for the massive case [SJ |H] apply to the massless case, and vice versa. 

It should be mentioned that, in the existing literature, the constant C used to be chosen as a 
real number different from zero. However, it is precisely this fact that, at the end, leads to the 
breakdown of the existing scalar reduction methods. Namely, for some kinematics (e.g. collinear 
on-shell external lines) the system IlKSt has no solution for C ^ 0. However, if the possibility C = 
is allowed, the system (|13|) will have a solution regardless of kinematics. This makes it possible 
to obtain additional reduction relations and formulate methods applicable to arbitrary number of 
external lines and to arbitrary kinematics. 

If (|12|l is taken into account, and after using the relation 

- Ej =1 ^o( D ^ {"* + Ski}) = (^y'loiD - 2; H}), (14) 
which can be easily proved from the representation (0), the relation (|1U|I reduces to 

N 

CI^{D-2-{ Vk \) = ^2 Zi I^(D-2;{u k -6 ki }) 

i=l 

N 

+ (^ 2 )(D - 1 - E »j) z oIo(D; K}), (15) 

where Zi are given by the solution of the system l|12f> . and the infinitesimal part proportional to ie has 
been omitted. This is a generalized form of the recursion relation which connects the scalar integrals 
in a different number of dimensions [3J0J|51|Sl[7j. The use of the relation (|15l) in practical calculations 
depends on the form of the solution of the system of equations <|12fl . For general considerations, it 
is advantageous to write the system (|12fl in the following way: 



( 


1 


1 


1 


\ 


/ 


-c \ 


( 


zo 


\ 


1 





ri2 ■ 


• r 1N 






Z\ 









1 


r\2 


• 


■ r 2 N 






Z2 









V i 


tin 


r 2 N ■ 


■ 


/ 


V 


ZN j 


V 





/ 



(16) 



G 



In writing we have taken into account the fact that zq = ~Y^ =l z^. In this way, the only free 
parameter is zq and by choosing it in a convenient way, one can always find the solution of the above 
system and, consequently, be able to use the recursion relations i|15f) . 

In the literature, for example in Ref. 00 [7], the recursion relations are obtained by inserting 
the general solution of the system l|12l) . i.e. the system (|16|1 . into the relation i|15|) . The recursion 
relations thus obtained are of limited practical use if the matrices of the mentioned systems are very 
singular. This happens when there arc cither two or more collinear external lines or, in general, 
for N > 6. When this is the case, the analysis of the general coefficient of the recursion becomes 
very complicated and in many cases unmanageable. There are cases when all coefficients vanish. As 
stated in for N > 7, owing to the drastic reduction of the recurrence relations these cases need 
a separate investigation. In addition, the above-mentioned problems with C ^ appear. To avoid 
these problems, a different approach to recursion relations can be taken. It is based on the fact, that 
finding any solution of the systems of equations mentioned above, makes it possible to perform the 
reduction. Being forced to use computers, it is very convenient and important that the reduction 
procedure be organized in a such a way that the recursion relations are classified and used depending 
on the form of the solutions of the above systems. If this is done, the increased singularity of the 
kinematic determinants turns out to be working in our favour by making it easy to find a solution 
of the systems of linear equations relevant to the reduction. 

In the following we frequently refer to two determinants, for which we introduce the notations: 
for the determinant of the system (|13fl we introduce det (Rn), while for the determinant of the system 
1 01) we use det(S , Ar). Depending on whether the kinematic determinants det(i?/v) and det(<S/v) are 
equal to zero or not, we distinguish four different types of recursion relations following from ||TSj. 
Before proceeding to consider various cases, note that in the case when det(S') ^ 0, it holds 

det(R N ) 

C -- Z °d^)- (17) 

It should be mentioned that for some of the recursion relations presented below one can find 
similar expressions in the literature. For reasons of clearness, connections of the relations given 
below with those existing in the literature are commented upon after the analysis of all possible 
cases has been considered. 

Let us now discuss all possible cases separately. 

4.1 Case I: det(S N ) ^ 0, det(R N ) ^ 

The most convenient choice in this case is Zq = 1. It follows from l|17l) that C ^ 0, so that the 
recursion relation 1|15|) can be written in the following form: 



4wn 2 (D- 1 -£ i=1 i/j) 



CI»(D-2;{v k }) 



V N zJ^{D-2-{v k -5 k% }) 

'l = 1 



(18) 



As it is seen, this recursion relation connects the scalar integral in D dimensions with the scalar 
integrals in D — 2 dimensions and can be used to reduce the dimensionality of the scalar integral. 

Since det (R n) ^ 0, some more recursion relations can be directly derived from (|1U|) . By directly 
choosing the constants Zi in (|10|) in such a way that Zi = for k = 1, • • • , N, we arrive at a system 
of N equations which is always valid: 

N N 

- rjfujl^D; {ui + %}) = ]T vjI?(D; + % - S lk }) 

N 

-(D-^i/XtAW), k = l,---,N. (19) 
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In the system (|19|l we have again disregarded the non-essential infinitesimal term proportional to 
ie. The matrix of the system (|19|l is the same as the matrix of the system <|12|) . whose determinant 
is different from zero, so that the system <|19|) can be solved with respect to I^(D;{vi + <%}), 
j = 1, ■■■,N. The solutions represent the recursion relations which can be used to reduce the 
powers of the propagators in the scalar integrals. Making use of these relations and the relation 
ijTSjl . each scalar integral 1^ (D; {i/i}) belonging to the type for which det(Sjv) ^ 0, det(J?jv) 
can be represented as a linear combination of integrals I^(D'; {1}) and integrals with the number of 
propagators which is less than N. For the dimension D', one usually chooses 4 + 2e, where e is the 
infinitesimal parameter regulating the divergences. Even in the case when one starts with D < D' , 
one can make use of the recursion (|18() to change from the dimension D to the dimension D' . 

In addition to the two sets of recursion relations presented above, by combining them one can 
obtain an additional and very useful set of recursion relations. This set at the same time reduces D 
and Vi in all terms. By adding and subtracting the expression Sj^Iq (D; {fj + 5ij — do-) in the first 
term on the right-hand side of the system H19|) and makeing use of the relation l|14|l . one finds 



N 



^(r fe -r J )Vo W (A{^ + <M) = -(4ir^r 1 l£ I (D-2;{v l -5 lk }) 

3 = 1 

N 

-(£-1-^)1^ (AM). k = l,.-.,N. (20) 
i=i 

The solution of this system of equations can in principle be used for reducing the dimension of the 
integral and the propagator powers. However, a much more useful set of the recursion relations is 
obtained by combining (|20[) and (JTSJ- Expressing the second term on the right— hand side of (|20|) 
with the help of l|18|) . leads to 

N r N 

^fa-r^Vo^S + = (4V)" 1 5>i - ^k)^ {D - 2; { Vi - ^}) 



3=1 

-CI»(D-2;{v i }) 



k = l,---,N, (21) 



where Z{ and C represent solutions of the system (fTB|) for zq = 1. Solutions of the system l|2"T|) rep- 
resent the recursion relations which, at the same time, reduce (make smaller) the dimension and the 
powers of the propagators in all terms (which is very important). As such, they are especialy conve- 
nient for making a rapid reduction of the scalar integrals which appear in the tensor decomposition 
of high-rank tensor integrals. 

4.2 Case II: det^) ^ 0, det(R N ) = 

The most convenient choice in this case is zq = 1. Unlike in the preceding case, it follows from (|17|) 
that C = 0, so that the recursion relation (|15|l can be written as 



4irfM 2 {D-l-J2 J= iV3) 



z i I»(D-2;{v k -5 ki }) 

* J l — l 



(22) 



It follows from l|22|l that it is possible to represent each integral of this type as a linear combination 
of scalar integrals with the number of propagators being less than TV. 

4.3 Case III: det(SW) = 0, det(R N ) ^ 

This possibility arises only if the first row of the matrix of the system l|16f) is a linear combination 
of the remaining rows. Then, the system l|ltj|) has a solution only for the choice zq = 0. With this 
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choice, the remaining system of equations reduces to the system <|12|) . where the constant C can be 
chosen at will. After the parameter C is chosen, the constants Zi are uniquely determined. Thus 
the recursion relation (|15|l with the choice C = 1 leads to 



I»(D; { Vk }) = V ' Zi I»(D; - 



(23) 



Consequently, as in the preceding case, the scalar integrals of the type considered can be represented 
as a linear combination of scalar integrals with a smaller number of propagators. 



4.4 Case IV: det(SV) = 0, det(R 







Unlike in the preceding cases, in this case two different recursion relations arise. To derive them, we 
proceed by subtracting the last, (N + 1)— th, equation of the system (|16fl from the second, third,... 
and N— th equation, respectively. As a result, we arrive at the following system of equations: 



/ o 





V i 



t 

ri2 - ri N 

n,N-i - nN 
tin 



1 

ru - r 2 N 
-r 2 N 



T2,N-l ~ r 2N 
T2N 



1 

tin 

T2N 



TN~l,N 





Zl 
'-2 



ZN-1 
Z N 



( *0 \ 






V o / 



(24) 



As it is seen, the first N equations of the above system form a system of equations in which the 
constant C does not appear, and which can be used to determine the constants Zi, i = 1, • • • N. 
The fact that det(Sjv) = implies that the determinant of this system vanishes. Therefore, for the 
system in question to be consistent (for the solution to exist), the choice zq = has to be made. 
Consequently, the solution of the system, Zi (i = 1, 2, ...-/V), will contain at least one free parameter. 
Inserting this solution into the last, (N + 1)— th, equation of the system (|24|l . we obtain 



r iN Zi = C . 



(25) 



Now, by arbitrarily choosing the parameter C, one of the free parameters on the left-hand side can 
be fixed. 

Sometimes, (for instance, when there are collinear external lines) the left-hand side of Eq. (|25|) 
vanishes explicitly, although the solution for zi contains free parameters. In this case the choice 
C = has to be made. 

Therefore, corresponding to the case when det(S , Ar)=det(i?Ar) = 0, one of the following two 
recursion relations holds: 



obtained from (|15|l by setting z$ = and C = 1, or 

° = E JV ' Vo W (Ah-4i}), 

z — ' I— 1 



Ski}), 



(26) 



(27) 



obtained from (|15|l by setting zo — and C = 0. 

In the case Ij26|l . it is clear that the integral with TV external lines can be represented in terms of 
the integrals with N — 1 external lines. What happens, however, in the case (|27(l ? With no loss of 
generality, we can take that Z\ ^ 0. The relation l|27|l can then be written in the form 



zJ^{D; { Vk }) = - zJ^(D; {v k + 5 kl - 5 kl }). 



(28) 
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We can see that, in this case too, the integral with N external lines can be represented in terms of 
the integrals with N — 1 external lines. In this reduction, v * remains conserved. 

Based on the above considerations, it is clear that in all the above cases with the exception of that 
when det(5jv) ^ 0, det(i?jv) ^ 0, the integrals with N external lines can be represented in terms of 
the integrals with smaller number of external lines. Consequently, then, there exists a fundemantal 
set of integrals in terms of which all integrals can be represented as a linear combination. 

Before moving on to determine a fundamental set of integrals, let us briefly comment on the 
recursion relations for scalar integrals that can be found in the literature. As we see below, det(>S/v) 
is proportional to the Gram determinant. All recursion relations for which the Gram determinant 
does not vanish are well known. Thus the relation of the type (|18|l can be found in Refs. 0^|S1|21|7], 
while the solutions of the systems (11911 and H21|) correspond to the recursion relations (28) and (30), 
respectively, given in Ref. [Hj • Even though Case II also belongs to the class of cases for which the 
Gram determinant is different from zero, the system (|13fl has no solution for C^O. This is a reason 
why the problem with using recursion relations appear in all approaches in which it is required that 
C^O. This can be seen from the discussion in [7J (the method is based on the choice C = 1) where 
the authors state that the reduction cannot be done for N — 3 with on-shcll external lines, and for 
N = 4 when one of the Mandelstam variables s or t vanishes. Such cases, however, are unavoidable 
when obtaining leading twist NLO PQCD predictions for exclusive processes at large momentum 
transfer. On the other hand, in the approach of Ref. UJ, where all coefficients of the recursion are 
given in terms of det(Sjv) and the minors of the matrix Sn, the relation of the type (|22|) can be 
obtained (Eq. (35) in Ref. |HJ). Cases III and IV, for which the Gram determinant vanishes are of 
special interest. One of the most discussed cases in the literature, belonging to Case III, is N = 6. 
The recursion relations of the type (|23[) can be found in Ref. ^ 03 El [7j. As for Case IV, it is 
especially interesting owing to the fact that it includes all cases for N > 7. In this case, the systems 
(|13fl and (|16fl have no unique solution. Case IV causes a lot of trouble for approaches in which the 
recursion coefficients are given in terms of det(Sjv) and the minors of the matrix Sn, for example 
in Ref. 0. The problem consists in the fact that all determinants vanish, making it impossible to 
formulate the recursion relation, so these cases need a separate investigation. On the other hand, 
the method of Ref. [7J, based on using pseudo-inverse matrices, can be used to construct the most 
general solution of the system l|13f) for the case of the vanishing Gram determinant. Even though 
the authors of Ref. claim that using their approach one can always perform the reduction of 
the V-point function (N > 6), that does not seem to be the case. Namely, the method in [7] is 
based on the choice C = 1, and as it has been shown above, in some cases belonging to Case IV 
the system (112(1 has no solution for C ^ 0, implying that the reduction cannot be performed. The 
impossibility of performing the reduction manifests itself such that v ■ K ■ v = (see (15) and (19) 
in [7]), a consequence of which is that the recursion coefficients become divergent. The situation 
of this kind arises regularly when dealing with integrals containing collincar external lines, i.e. for 
exceptional kinematics. The method of Ref. [7J has been obtained for non-exceptional kinematics. 

In view of what has been said above, most of the problems with existing reduction methods 
appear when dealing with the integrals with a large number of external lines. In all considerations 
in the literature that happens for N > 6. It is very important to point out that, this is valid for 
the case when external momenta span the 4-dimensional Minkowski space. If the dimensionality 
of space span by external momenta is smaller, the problems start appearing for smaller N. Even 
though one can find the statements that such cases are at the moment of minor physical interest, 
we disagree. Namely, as stated earlier, the analysis of exclusive processes in PQCD, even for simple 
processes, requires evaluation of the diagrams with N > 6. Since these diagrams contain collinear 
external lines, the kinematics is limited to (d < 4)-dimensional subspace. Thus, for example, for 
nucleon Compton scattering the integrals with N — 8 external lines contribute and the kinematics is 
limited to the 3-dimensional subspace. A consequence of this is that the problem with using existing 
reduction methods will start appearing at the level of one-loop N — 5 diagrams. 

The reduction method presented in this paper is formulated with an eye on exclusive processes 
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in PQCD. The main point of the method is that the reduction is defined in terms of the solution of 
the linear systems given by (|13fl and (|16fl . A consequence of this is that the method is quite general, 
very flexible, practical and easily transfered to the computer program. To perform reduction, one 
only needs to find solution of the above systems which can always be done. A very pleasing feature 
of this reduction is that the increased singularity of the kinematic determinants facilitates reduction, 
since finding a solution of the relevant linear systems becomes easy. 

5 On the fundamental set of integrals 

We now turn to determine the fundamental set of integrals. To this end, let us first evaluate the 
determinant of the system Ijl6|l . det(5W), and determine the conditions under which this kinematic 
determinant vanishes. 

By subtracting the last column from the second, third, ... and TV— th column, respectively, and 
then the last row from the second, third, ... and N— th row, respectively, we find that the determinant 
det(S'Ar) is given by the following expression: 

det(S N ) = -det[-2(n - r N )(rj -r N )}, i,j = 1, . . . , N - 1. (29) 

As it can be seen, det(,5/v) is proportional to the Gram determinant. Denote by n the dimension of 
the vector space spanned by the vectors rj — rjy, (i — 1, . . . , N — 1). Owing to the linear dependence 
of these vectors, the determinant vanishes when N > n + 1. As in practice, we deal with the 
4— dimensional Minkowski space, the maximum value for n equals 4. An immediate consequence of 
this is that all integrals with N > 5 can be reduced to the integrals with N < 5. 

In view of what has been said above, all one-loop integrals are expressible in terms of the integrals 
1^(4 + 2e;{l}), j£(4 + 2e;{l}), 7^(4 + 2e;{l}), belonging to Case I, and the general two-point 
integrals Iq(D' ; v\,V2)i which are simple enough to be evaluated analytically. 

Next, by substituting D = 6 + 2e, N = 5 and z/j = 1 into the recursion relation Q15p. one finds 

5 

C/ 5 (4 + 2e; {1}) = ^ zJ^A + 2e; {S kk - S ki }) + (47r M 2 )(2e) z Q I 5 (6 + 2e; {1}). (30) 

i=l 

Owing to the fact that the integral Iq (6 + 2e; {1}) is IR finite^], the relation (|3L)fl implies that the 
TV = 5 scalar integral, Iq (4 + 2e; {1}), can be expressed as a linear combination of the N = 4 scalar 
integrals, Iq (4 + 2e; {1}), plus a term linear in e. In massless scalar theories, the term linear in e 
can simply be omitted, with a consequence that the TV = 5 integrals can be reduced to the N = 4 
integrals. On the other hand, when calculating in renormalizable gauge theories (like QCD), the 
situation is not so simple, owing to the fact that the rank P(< iV) tensor integrals are required. 

In the process of the tensor decomposition and then reduction of scalar integrals all way down 
to the fundamental set of integrals, there appears a term of the form (1/£)/q(4 + 2e;{l}), which 
implies that one would need to know an analytical expression for the integral Iq(4 + 2e;{l}), to 
order e. Going back to the expression (|30|> . we notice that all such terms can be written as a 
linear combination of the box (4— point) integrals in 4 + 2s dimensions and 5— point integrals in 

6 + 2e dimensions. Therefore, at this point, the problem has been reduced to calculating the integral 
Iq(6 + 2e;{l}), which is IR finite, and need to be calculated to order O(e ). It is an empirical 
fact 0] 03 El El that in final expressions for physical quantities all terms containing the integral 
ig(6 + 2e;{l}) always combine so that this integral ends up being multiplied by the coefficients 
0(e), and as such, can be omitted in one— loop calculations. A few theoretical proofs of this fact 
can be found in literature 0JEOE!j but, to the best of our knowledge, the proof for the case of 
exceptional kinematics is still missing. That being the case, in concrete calculations, (to be sure 
and to have all the steps of the calculation under control ), it is absolutely necessary to keep track 
of all the terms containing the integral Jq (6 + 2e; {!}), add them up and check whether the factor 
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multiplying it is of order O(e). Even though the experience gained in numerous calculations shows 
that this is so, a situation in which the integral Iq (6 + 2e; {1}) would appear in the final result for 
a physical quantity accompanied by a factor 0(1) would not, from practical point of view, present 
any problem. Namely, being IR finite, although extremely complicated to be evaluated analyticaly, 
the integral Iq (6 + 2e; {1}) can always, if necessary, be evaluated numerically. 

Based on the above considerations we may conclude that all one— loop integrals occurring when 
evaluating physical processes in massless field theories can be expressed in terms of the integrals 

I*{D';v x ,v 2 ), 7 3 (4 + 2 £ ; 1,1,1), I 4 (4 + 2e; 1, 1, 1, 1). 

These integrals, therefore, constitute a minimal set of fundamental integrals. 

In view of the above discussion, we conclude that the set of fundamental integrals is comprised of 
integrals with two, three and four external lines. Integrals with two external lines can be calculated 
analytically in arbitrary number of dimensions and with arbitrary powers of the propagators. They 
do not constitute a problem. As far as the integrals with three and four external lines are concerned, 
depending on how many kinematic variables vanish, we distinguish several different cases. We now 
show that in the case N = 3 we have only one fundamental integral, while in the case corresponding 
to N = 4 there are six integrals. For this purpose, we make use of the vanishing of the kinematic 
determinants det(i?jv) and det(Sjv). 



Jg(J;^)^rN ' ■ (3D 



5.1 The general scalar integral for N = 2 

According to iJSJ, the general massless scalar two— point integral in D space— time dimensions is of 
the form 

d D l 1_ 

(2tt) d A? A? 

The closed form expression for the above integral, valid for arbitrary D — n + 2s, and arbitrary 
propagator powers v\ and v 2 , is given by 

I%{n + 2e; v x , v 2 ) = (4^ 2 ) 2 -/ 2 (-l)-+- (-p 2 - ie )"/ 2 -— 

T (y x +v 2 - n/2 - e) T (n/2 -v l +e)T {n/2 - v 2 + e) 

x r(-e) r(i + e) r(i + E ) 

/ 2 (4 + 2e,l,l), (32) 



1 T(2 + 2 £ ) 



T{v{)T{y2)T{n-vi-V2 + lE) 
where 

i, ^ if p 2 +ie\ e r(-e)r 2 (l+e) , . 

'o^M) = (4^(-^j l r(2 + 2 £ ) " ^ 

It is easily seen that in the formalism of the dimensional regularization the above integral vanishes 
for p 2 = 0. 



5.2 The scalar integrals for N = 3 

The massless scalar one— loop triangle integral in D = 4 + 2e dimensions is given by 

7 »' 4 + 2e '( 1 »-<" 2 ^*/(S^3^' ,34) 

Making use of the representation JJJ, and introducing the external masses pf = mf (i = 1,2,3), the 
integral (|34(l can be written in the form 



/ dxidx2da;3 5(xi + x 2 + x 3 — 1) 
Jo 



x (—x\X2 m\ — X2X3 to 2 — X3X1 ml — ie) e . (35) 
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It is evident that the above integral is invariant under permutations of external massess mf. De- 
pending on the number of the external massless lines, and using the above mentioned symmetry, 
there are three relevant special cases of the above integral. We denote them by 

I 1 ™ = 7 3 (4 + 2e,{l}; 0,0,7^), (36) 
I 2m ee 7 3 (4 + 2e,{l};0,m2,m2), (37) 
Il m ee ll{A,{l}-mlmlml), (38) 

The integrals I 3 m and 7f m are IR divergent and need to be evaluated with e > 0, while the integral 
if™ is finite and can be calculated with e = 0. 

Now, it is easily found that the determinants of the systems of equations l|13|) and (|16|l are, for 
N = 3, given by 

det(i? 3 ) = 2m\m\m\, (39) 
det(S 3 ) = (mf) 2 + (m\) 2 + (ml) 2 - 2m\m\ - 2m\m\ - 2m\m\, (40) 

As is seen from (|39ll . if at least one of the external lines is on mass-shell, the determinant det(i? 3 ) 
vanishes. Consequently, using the recursion relations (Case II or IV) the integrals l\ m and I 2m can 
be reduced to the integrals with two external lines. Therefore, we conclude that among the scalar 
integrals with three external lines the integral if™ 1 is the only fundamental one. 

The result for this integral is well known ^1 ^3 ^] . In it is expressed in terms of the 
dimensionless quantities of the form 



1 

x i.2 = ^ 




l--4 + -#±W l--4-^-4 -4 



m 2 



(41) 



and, being proportional to l/(x\ — x 2 ), appears to have a pole at x\ = x 2 . It appears that the final 
expression |12j is not well defined when x\ = x 2 . 
On the basis of Eas. ll4T)|l and ll4*T|l . one finds that 



1 



Xi — X2 = —y/det(S 3 ). 

This equation implies that when x\ — X2 = 0, instead of examining the limit of the general expression 
in one can utilize the reduction relations 11' Ml) (corresponding to det(i? 3 ) 7^ and det(S' 3 ) = 0) 
to reduce the IR finite integral with three external lines, if" 1 , to the integrals with two external 
lines. 

5.3 The scalar integrals for N = A 

The massless scalar one— loop box integral in D = 4 + 2e space— time dimensions is given by 

<(« + »,{!}) -(MV I ^-AaUa, - < 42 » 

Making use of (|7|), introducing the external "masses" p 2 = m\ (i = 1,2,3,4), and the Mandelstam 
variables s = (p\ + p 2 ) 2 and t = (p 2 +p 3 ) 2 , the integral (j4*2")l becomes 

i r(2-e) f 1 



Jp (4 + 2s, {1}) — — ^— / dxidx 2 dx 3 dx4 S(xi + x 2 + x 3 + x 4 - 1) 

(Att) z (Att^Y J 



(— x\x 3 t — x 2 xa s — x\x 2 m\ ~x 2 x 3 m 2 — x 3 xa m\ — x\X\ m\ — ie) £ 2 . (43) 
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Introducing the following set of ordered pairs 

(s,t), (m 2 ,m 2 ), (m\,m\), (44) 

one can easily see that the integral (|43|l is invariant under the permutations of ordered pairs, as well 
as under the simultaneous exchange of places of elements in any two pairs. 

The determinants of the coefficient matrices of the systems of equations (|13|1 and (|16|l . corre- 
sponding to the above integral, are 

det(i? 4 ) = s 2 t 2 + {mjml) 2 + (m 2 m 2 ^) 2 

— 2stm\m\ — 2stm\m\ — 2m\m\rn\m\ . (45) 



det(S 4 ) = 2 [ st (ml 



+ m\m\(s + t + ml + m 



i ™2„2/ i , 2 

+ m 1 m 3 (s + t — m 1 - 
— s(mlm\ + m^ml) 



ml 

2 
3 

2 , 2 



m, 



t) 

-ml) 



3 T ni-2 T W 4 ) 
j.1 2 2 i 2 2\ - 

t(m 1 m^ + m 2 m 3 ) 



(46) 



By looking at the expression for det(i? 4 ) given in (|45|l it follows that all box integrals Iq that are 
characterized by the fact that in each of the ordered pairs in (I44|l at least one kinematic variable 
vanishes, are reducible. Therefore, for a box integral to be irreducible, it is necessary that both 
kinematic variables in at least one of the ordered pairs should be different from zero. Owing to the 
symmetries valid for the box integrals it is always possible to choose that pair to be (s,t). 

Taking into account symmetries, and the number of external massless lines, there are six poten- 
tially irreducible special cases of the integral Adopting the notation of Ref. 0], we denote 
them by 



j4m 
1 i 
j3m 
1 4 

r2mh 
j2me 



rim, 
1 4 
rOrn 
1 4 



J o (4, {1}; 8, t, ml, ml, m%, ml), 
7 4 (4 + 2s, {1}; s, t, 0, m\, m\, m 2 ), 
I*(4 + 2£,{l};s,t,0,0,m 2 3 ,m 2 ), 
l£(4 + 2e,{l}; S ,t,0,m 2 2 ,0,ml), 
l£(4 + 2e,{l}; S ,t,0,0,0,m 2 ), 
7 4 (4 + 2 £ ,{l}; S ,i,0,0,0,0,), 



(47) 
(48) 
(49) 
(50) 
(51) 
(52) 



with all kinematic variables appearing above being different from zero. The results for these integrals 
are well known gl Q31 [TJ Ej . 

The integrals Q48I52H are IR divergent, and as such need to be evaluated with e > 0, while the 
integral l|47|l is finite and can be calculated in D = 4. The results for these integrals, obtained 
in ^2 EH f° r arbitrary values of the relevant kinematic variables, and presented in a simple and 
compact form, have the following structure: 



IZ(s,t;mi) 



r(l-e)r 2 (l + e) 



(4tt) 2 T(l + 2e) 
r G A '(s,t;e;m 2 ) 



x/dit(lf) 
+ H K ( S ,t;m 2 ) 



0(e), 



K € {0m, lm, 2me, 2mh, 3m, 4m} . 



(53) 



The IR divergences (both soft and collinear) of the integrals are contained in the first term within 
the square brackets, while the second term is finite. The function G K (s,t; e;mf) is represented by 
a sum of powerlike terms, it depends on e and is finite in the e — > limit. As for the function 
H K (s,t;m 2 ), it is given in terms of dilogarithm functions and constants. In the above, det(Rf) is 
the determinant corresponding to the integral I± given in (|4 71521) . 
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For the purpose of numerical integration, it is very useful to have the exact limit of the integral 
Ijf when det(i?4") — * 0. This limit can be determined in an elegant manner by noticing that for 
det(Rf ) — the reduction relations corresponding to Cases II and IV apply, making it possible 
to represent the box integral as a linear combination of the triangle integrals. This result can 
be made use of to combine box and triangle integrals (or pieces of these integrals) with the aim of 
obtaining numerical stability of the integrand [2]. 

The integrals (|47I52|) are irreducible only if the corresponding kinematic determinant det(R^) 
does not vanish. 

With the help of the tensor decomposition and the scalar reduction procedures, any dimensionally 
regulated one— loop N— point Feynman integral can be represented as a linear combination of the 
integrals: 

I%(D';v u v 2 ), 

r3m 
J 3 i 

r4m r3m jlmh jlrae j\ra jQm (KA\ 

4 ' 4 5 4 j J 4 i J 4 i 1 4 i \ oq: ) 

multiplied by tensor structures made from the external momenta and the metric tensor. The inte- 
grals in (|54|) constitute a fundamental set of integrals. An alternative and more convenient set of 
fundamental integrals is obtained by noticing that all the relevant box integrals are finite in D = 6. 
On the basis of Eq. Ijl5|l . all IR divergent box integrals can be expressed as linear combinations 
of triangle integrals in D = 4 + 2e dimension and a box integral in D = 6 + 2e dimension. Next, 
using the same equation, all triangle integrals can be decomposed into a finite triangle integral and 
two— point integrals. In the final expression thus obtained all divergences, IR as well as UV, are 
contained in the general two-point integrals and associated coefficients. Therefore, an alternative 
fundamental set of integrals is comprised of 

lt(D'- Vl ,v 2 ), 

r3m 
i 3 J 

r4m j3m j2mh j2me jlm jOm fKK\ 

4 ' 4 5 4 ' A ' A ! 4 ' \ 00 ) 

where denotes box integrals in D — 6 dimensions, explicit expressions for which are given in the 
Appendix A. A characteristic feature of this fundamental set of integrals, which makes it particularly 
interesting, is that the integral Iq is the only divergent one, while the rest of integrals are finite. 



6 Conclusion 

In this work we have considered one— loop scalar and tensor Feynman integrals with an arbitrary 
number of external lines which are relevant for construction of multi— parton one— loop amplitudes 
in massless field theories. 

Main result of this paper is a scalar reduction approach by which an arbitrary A— point scalar 
one— loop integral can be reqursively represented as a linear combination of eight basic scalar integrals 
with rational coefficients depending on the external momenta and the dimensionality of space— time, 
provided the external momenta are kept in four dimensions. The problem of vanishing of the 
kinematic determinants, which is a reflection of very complex singularity structure of these integrals, 
has been solved in an elegant and transparent manner. Namely, the approach has been taken 
according to which instead of solving the general system of linear equations given in (|12J) . and then 
finding the limit, which sometimes doesn't exists, of the obtained solution corresponding to a given 
singular kinematic situation, we first obtain and then solve the system of equations appropriate to 
the situation being considered. 

Our method has been derived without any restrictions regarding the external momenta. As 
such, it is completely general and applicable for arbitrary kinematics. In particular, it applies 
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to the integrals in which the set of external momenta contains subsets comprised of two or more 
collincar momenta. This kind of integrals are encountered when performing leading— twist NLO 
PQCD analysis of the hadronic exclusive processes at large— momentum— transfer. Trough the tensor 
decomposition and scalar reduction presented, any massless one-loop Feynman integral with generic 
4-dimensional momenta can be expressed as a linear combination of a fundamental set of scalar 
integrals: six box integrals in D — 6, a triangle integral in D — 4, and a general two— point integral. 
All the divergences present in the original integral are contained in the general two-point integral 
and associated coefficients. 

In conclusion, the computation of IR divergent one— loop integrals for arbitrary number of ex- 
ternal lines can be mastered with the reduction formulas presented above. The iterative structure 
makes it easy to implement the formalism in algebraic computer program. With this work all the 
conceptual problems concerning the construction of multi— parton one— loop amplitudes are thus 
solved. 
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Appendix A 

In addition to the explicit calculation, the irreducible box integrals in D = 6 dimensions can be 
obtained using the existing analytical expressions for the irreducible box integrals in D = 4 + 2e 
dimensions and the reduction formula (|15|) . To this end, we substitute D = 6 + 2s, N = 4, Vi = 1 
and C = 1 into the relation l|15(l and find 

7 4 (6 + 2s; {1}) = +1) Zq Ko(4 + 2e; {1}) - £ z 4 / 4 (4 + 2e; {S kk - 6 ki })^j . (56) 

Note that the IR divergences in D = 4 + 2e box integrals are exactly cancelled by the divergences 
of the triangle integrals. 

The expressions for the relevant basic massless scalar box integrals in D — 6 space— time dimen- 
sions are listed below: 
The three— mass scalar box integral 




2 



t + ie 



t — m 2 — m 3 + 2m\m\ 



1 - (t + ie) f m 
1 - K + ie) f m 

I 3 (ml,ml,t) 



t - m? 



st — m^m 2 i 



In 



t + ie 
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— — ( s — m\ — m\ + 2m\m\- 



l st — m\m\ 

The adjacent ("hard") two— mass scalar box integral 



T 3 (7713,7714, s) 



ll mh (D = 6;s,t;mi,mi) = 



i 1 



(47r) 2 47T/X 2 



7773 + ie 



777^ + it 



Li 2 1 



7774 + ie 
s + ie 



Ll 2 1 - T~ . + Ll 2 



+ Li 2 



t + ie 
1 - (t + ie) / 



2/7/7; 



- Li 2 

2 2 \ 
77737774 \ 



- ^s- 777^ -ml + 2 ; , 



1 - (s + ie) / 2mfc 
1 - (ml + ie) f 2mh 

I 3 (ml,ml,s) }. 



The opposite ("easy") two— mass scalar box integral 



J 4 2me (73 = 6;s,t; m 2 ,™ 2 ) 



i 1 



(47r) 2 47T/Z 2 



x /7 2me <^ Li 2 



1 - (s + ie) / 



2m e 



Li? 



l-(t + ie)/ 2 



- Li 2 



1 - (m 2 + ie) / 



2me 



- Li2 



1 - (777 2 + ie) / 



2me 



The one— mass scalar box integral 



ll m (D = 6;s,t;ml) = 



(47r) 2 47T/J 2 



x h lm { Li 2 



l-^ + ie)/ 1 



+ Li 2 



l-^ + ie)/ 1 



Li 2 



1 - (ml + ie) / 



lm 



7T 



The zero— mass (massless) scalar box integral 



I» m (D = 6;s,t) 



x h° m l Li 2 



(4vT) 2 47T/J 2 

l-(.s + ie)/ 



Dm 



+ Li 2 



1 - (t + ie) /' 



Ora 



IT 



where 



^ = ( oVdeW) 



det(5 4 K ) 
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and 



(4tt) 2 

The functions appearing above are given by 



X 3 (o, b, c) = I$ m (D = 4; a, b, c) . (63) 



j3m j2me 

jlrah ^lm 



S + t — mf, — m| 

st — m\m\ 
s + t — m| 

s + t 



st 



,3m ( , . 2 2 2 . 2 mjt + mjs - 2mjm 4 

a = s + t — m 9 — rrir, — m A + m-, 5 — 5 

V 2 01 4 6 st - ra\m\ 

/ 2 2 

/i 2me = (s + t-mi-mty\ 
h lm = (s + t-miy 1 , 
h° m = {s + ty 1 . 

Appendix B 

As an illustration of the tensor decomposition and scalar reduction methods, we evaluate an one-loop 
6-point Feynman diagram shown in Fig. [3 

Note that, due to the kinematics which is bounded to three dimensional Minkowski subspace 
there are no four linearly independent four- vectors. Consequently, this diagram is of complexity of 
the 7-point one-loop diagram with four dimensional external kinematics. We choose this particular 
diagram because of the compactness of intermediate and final expressions. 

This is one (out of 462) diagram contributing to the NLO hard-scattering amplitude for the exclu- 
sive process j(k\, £1) 7(^2, £2) —> 7T+ (P+) ^(P-), (with both photons on-shell) at large momentum 
transfer. 

In the 77 centre-of-mass frame, the 4-momenta of the incoming and outgoing particles are 

fci,2 = V5/2(l ) Tsin0 c .m.,O ) ±cos0 c . m .), P± = yfs/2 (1,0,0, ±1), (64) 

while the polarization states of the photons are 

ef =e 2 F =Tl/V2{0,cose c . m .,±i,sm9 c . m .), (65) 

where y/s is the total centre-of-mass energy of the 77 system (or the invariant mass of the tt + tt~ 
pair) . 

For example, taking Q c . m , — n/2 and assuming that the photons have opposite helicities, the 
amplitude corresponding to the Feynman diagram of Fig. [21 is proportional to the integral 

2 J I2 { l +pi) 2 {l „ pi) 2 {l+pb) 2 {l+pfi) 2 {l+p7) 2 > W 

with the momenta pi (i = 1, . . . , 7) 

Pi=xP + , p 2 = xP + , p 3 = fci-yP_, p A = yP_, 

Pb = yP-, Pe = yP--k 2 , P7 = xP+ + yP- - fa- (67) 
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Figure 2: One of the diagrams contributing to the hard scattering amplitude of the process 77 
7T+ tt- at NLO. 



The quantities x and x = 1 — m (y and y = 1 — y) are the fractions of the momentum P + (P_ ) shared 
between the quark and the antiquark in the 7r + (7r - )- 

With the aim of regularizing the IR divergences, the dimension of the integral is taken to be 
D = A + 2e. 

The integrall is composed of one-loop 6-point tensor integrals of rank 0, 1, 2, 3 and 4. Performing 
the tensor decomposition and evaluating the trace, we obtain the integral X in the form 

T = -2(1 + e) 2 [24s 3 xxyi/(47iV) 4 / ( ?(12 + 2e, {1,1, 1,1,1,5}) 



+ 6 s 3 y y (4 tt// 2 ) 7 6 (12 + 2 e, {1, 4, 1, 1, 2, 1}) 



s 3 y? 

+ 2 s 3 y (y - y) (4 tt ^ 2 )" / 6 (12 + 2 e, {1, 3, 2, 1, 2, 1}) 
+ 8 s 3 y y (4 tt/?) 4 7 6 (12 + 2 e, {1, 3, 1, 1, 3, 1}) 
+ 2s 3 xyy (4^ 2 ) 4 7 6 (12 + 2 e, {1, 2, 2, 2, 1, 2}) 
+ s 3 y (y-y) (4 7r M 2 ) 3 / 6 (10 + 2e,{l,2,2,l,2,l}) 
+ s 2 yy (1 + e) (4tt M 2 ) 7 6 (6 + 2e, {1, 1, 1, 1, 1, 1}) 
+ s (1 + e) (2 + e) (4tt/x 2 ) 2 ij(8 + 2e, {1,1, 1,1, 1,1}) 
+ ... 75 similar terms 



(68) 



Next, performing the scalar reduction using the method described in the paper, we arrive at the 
following expression for the integral written in terms of the basic integrals: 



1 = 



8 (1 + ef { (4 vr fi 2 ) [J 7 4 lm (6 + 2 e; -a/2, -s y/2; -s y/2) 
+ i±£ ji-( 6 + 2 e; -s y/2, -s/2; -a y/2) 

+ (l + e(l-|)) ll m (6 + 2e--sx/2,-sy/2;-s(xy + xy)/2) 



+ 



V a; 



el 



/ 4 2me (6 + 2 e; -saf/2, -s y/2; -a/2, -a (xy + x y)/2) 



f-T^-T + 2 - £ ] /2U + 2 : -.sT/2) 
(x — x)y \e (x — x) x, 



+ ( + 2 - % ] A, ( -I + 2 : - , , 2 1 

[x — x)y \e (x — x) x, 
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(y-y)x \e(y-y) 



2y 



2- 



I 2 (4 + 2e;-sy/2) 



2y — + 2-l) I 2 (4 + 2e;-sy/2) 



(y-y)x \e(y-y) y 

(1 — xy — 3yx)(l — yx — 3xy) 2(xy + xy)(8xxyy — xx — yy) 



xxyy(x- x)(y - y) 



exxyy(x - x) 2 (y - y) 2 



xl 2 (4 + 2e:~s(xy + x y)/2) + ixy + x _ y ^ h (4 + 2 e: -a/2) 

xxyy 



(69) 



Here, I 2 is the two-point scalar integral in D = 4 + 2e with i/$ = 1, while l\ m and /I™ 16 are box 
scalar integrals in D = 6 + 2e. Analytic expressions for these integrals are given in the Appendix 
A. Expanding Eq. (|69[) up to order O(e ), we finally get 



, , 1 -— Li 2 (a; - x) - ^Li 2 (a; - x) - -^ = Li 2 (y - y) - "-^Li 2 (y- y) 
(4 7r)^ s ( xy yx yx xy 

, xy + yx ir 2 xy + yx xy + xy ( s 
+ — =-=-Li 2 {-{x - x)(y - y)) + ^^ln — — 

xxyy o xxyy xxyy \2y 



(x-2x) 
x y (x — x) 
(V~2y) 

xy(y-y) 



In 



In 



In 2 



s x 
2~J^ 

sy 

2fi 2 

S X 

2^2 



y (x — x) 2 



(UL 

x(y-y) 2 *~ \2/j, 2 



(x-2x) 
xy (x — x) 

(y-2y) 
xy (y - y) 

X 



y (x — x) 2 

V 

x (y - y) 2 



La- 



in 
In 
In 2 



s x 
2^? 
s y 
2^? 
sx 
2~J? 



In* 



i[ sy 
2 y 2 



(1 — xy — 3yx)(l — yx — 3xy) f s (xy + xy) 



xxyy (x -x)(y -y) 
(xy + xy)(8xxyy - xx - yy) 2 
xxyy(x-x) 2 (y-y) 2 



2fi 2 
s (xy + xy) 



+ 



y (x — x) 2 

y 



In 



2^ 



x(y-y) 2 



In 



sy 
2y 2 



y (x — x) 2 
V 



In 



2^ 2 



2^ 



x(y- yf 



In 



s y 
~2yJ 



(xy + xy)(8xxyy — xx — yy) fs (xy + xy) 



xxyy(x - x) 2 (y - y) 2 
where 1/i = 1/e + 7 — ln(47r). 



2^2 



(70) 



References 

[1] G. Passarino and M. Veltman, Nucl. Phys. B 160 (1979) 151; G. J. van Oldenborgh and J. A. 
M. Vermaseren, Z. Phys. C 46 (1990) 425; W. L. van Neerven and J. A. M. Vermaseren, Phys. 
Lett. B 137 (1984) 241. 

[2] A. I. Davydychev, Phys. Lett. B 263 (1991) 107. 

[3] Z. Bern, L. J. Dixon and D. A. Ko sower, Phys. Lett. B 302 (1993) 299 [Erratum-ibid. B 318 
(1993) 649] |arXiv:hep-ph/9212308| . 



20 



Z. Bern, L. J. Dixon and D. A. Kosower, Nucl. Phys. B 412 (1994) 751 arXiv:hep-ph/9306240 
O. V. Tarasov, Phys. Rev. D 54 (1996) 6479 j arXiv:hep-th/9606018| . 

J. Fleischer, F. Jege rlehncr and O. V. Tarasov, Nucl. Phys. B 566 (2000) 423 
arXiv:hep-ph/9907327 . 

T. Binoth, J. P. Guillet and G. Heinrich, Nucl. Phys. B 572 (2000) 361 arXiv:hep-ph/9911342 ; 
G. Heinrich and T. Binoth, Nucl. Phys. Proc. Suppl. 89 (2000) 246 |arXiv:h ep-ph/0005324 

T. Binoth, J. P. Gui llet, G. Heinrich and C. Schubert, Nucl. Phys. B 615 (2001) 385 
|arXiv:hep-ph/0106243| . 

J. M. Campbell, E. W. Glover and D. J. Miller, Nucl. Phys. B 498 (1997) 397 
arXiv:hep-ph/9612413 . 



A. Denner and S. Dittmaier, Nucl. Phys. B 658 (2003) 175 arXiv:hep-ph/0212259 . 
G. Duplancic and B. Nizic, Eur. Phys. J. C 20 (2001) 357 arXiv:hep-p"h /0006249 . 
G. Duplancic and B. Nizic, Eur. Phys. J. C 24 (2002) 385 arXiv:hep-ph /0201306 . 
A. I. Davydychev, J. Phys. A 25 (1992) 5587. 

K. G. Chetyrkin and F. V. Tkachov, Nucl. Phys. B 192 (1981) 159; F. V. Tkachov, Phys. Lett. 
B 100 (1981) 65. 

F. Jegerlehner and O. Tarasov, Nucl. Phys. Proc. Suppl. 116 (2003) 83 arXiv:hep pli/02 12001 . 

A. T. Suzuki, E. S . Santos and A. G. Schmidt, Eur. Phys. J. C 26 (2002) 125 
|arXiv:hep-th/0205i58| . 

[17] C. Anastasiou, E. W . Glover and C. Oleari, Nucl. Phys. B 565 (2000) 445 
|arXiv:hep-ph/9907523j . A. T. Suzuki, E. S. Santos and A. G. Schmidt, |arXiv:hep-ph /0210083 



21 



